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Abstract 

. . The new release of data from Wilkinson Microwave Anisotropy Probe improves the observational status 

^ '• of relic gravitational waves. The 7-year results enhance the indications of relic gravitational waves in the 

{^JQI existing data and change to the better the prospects of confident detection of relic gravitational waves by the 

currently operating Planck satellite. We apply to WMAP7 data the same methods of analysis that we used 

V; . earlier [W. Zhao, D. Baskaran, and L.P. Grishchuk, Phys. Rev. D 80, 083005 (2009)] with WMAP5 data. 

'^ ' of consecutive WMAP data releases that the maximum likelihood value of the quadrupole ratio R, which 



in 



X 



characterizes the amount of relic gravitational waves, increases up to i? = 0.264, and the interval separating 



^^ . this value from the point i? = (the hypothesis of no gravitational waves) increases up to a 2cr level. The 

primordial spectra of density perturbations and gravitational waves remain blue in the relevant interval of 
wavelengths, but the spectral indices increase up to ng — 1.111 and nt = 0.111. Assuming that the maximum 
likelihood estimates of the perturbation parameters that we found from WMAP7 data are the true values 
of the parameters, we find that the signal-to-noise ratio S/N for the detection of relic gravitational waves 
by the Planck experiment increases up to S/N ~ 4.04, even under pessimistic assumptions with regard to 
residual foreground contamination and instrumental noises. We comment on theoretical frameworks that, 
in the case of success, will be accepted or decisively rejected by the Planck observations. 
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I. INTRODUCTION 

The Wilkinson Microwave Anisotropy Probe (WMAP) Collaboration has released the results of 
7-year (WMAP7) observations [l|, I4]. In this paper, we apply to WMAP7 data the same methods 
of analysis that we have used before 3] in the analysis of WMAP5 data. This is important for 
updating the present observational status of relic gravitational waves and for making more accurate 
forecasts for the currently operating Planck mission [4j]. 

In Sec. In] we briefly summarize our basic theoretical foundations and working tools. Part of 
this material was present in the previous publication [3|, but in order to make the paper self- 
contained we briefly repeat it here. Section IIIII exposes full details of our maximum likelihood 
analysis of WMAP7 data. In the focus of attention is the interval of multipoles 2 < i < 100, where 
gravitational waves compete with density perturbations. We compare all the results that we derived 
by exactly the same method from WMAP7, WMAP5 and WMAP3 datasets. This comparison 
demonstrates the stability of data and data analysis. On the grounds of this comparison, one can 
say that the perturbation parameters found from the consecutive WMAP data releases have the 
tendency of saturating near some particular values. The WMAP7 maximum likelihood (ML) value 
of the quadrupole ratio R is close to previous evaluations of R, but increases up to i? = 0.264. 
The interval separating this ML value from the point R = (the hypothesis of no gravitational 
waves) increases up to a 2a level. The primordial spectra remain blue, but the spectral indices in 
the relevant interval of wavelengths increase up to Ug = 1.111 and rit = 0.111. 

In Sec. IIVI we analyze why, to what extent, and in what sense our conclusions with respect 
to relic gravitational waves differ from those reached by the WMAP Collaboration. The WMAP 
team has found "no evidence for tensor modes." A particularly important issue, which we discuss 
in some detail, is the presumed constancy (or simple running) of spectral indices. We derive an 
exact formula for the spectral index rit as a function of wavenumbers and discuss in this context the 
formula for running that was used in WMAP analysis. Another contributing factor to the difference 
of conclusions is the difference in our treatments of the inflationary "consistency relations" based 
on the inflationary "classic result." We do not use the inflationary theory. 

A comprehensive forecast for Planck findings in the area of relic gravitational waves is pre- 
sented in Sec. |Vl We discuss the efficiency of various information channels, i.e. various correlation 
functions and their combinations. We perform multipole decomposition of the calculated S/N 
and discuss physical implications of the detection in various intervals of multipole moments. We 
stress again that the B-mode detection provides the most of S/N only in the conditions of very 



deep cleaning of foregrounds and relatively small values of R. The improvements arising from a 
28-month, instead of a nominal 14-month, Planck survey are also discussed. In the center of our 
attention is the model with the WMAP7 maximum likelihood set of parameters. For this model, 
the signal-to-noise ratio S/N in the detection of relic gravitational waves by Planck experiment 
increases up to S/N = 4.04, even under pessimistic assumptions with regard to residual foreground 
contamination and instrumental noises. Section IVII gives Bayesian comparison of different theo- 
retical frameworks and identifies predictions of R that may be decisively rejected by the Planck 
observations. 

II. PERTURBATION PARAMETERS AND CMB POWER SPECTRA 

The temperature and polarization anisotropics of CMB are produced by density perturbations, 
rotational perturbations and gravitational waves. Rotational perturbations are expected to be 
very small and are usually ignored, and are in this paper, too. The cosmological perturbations are 
characterized by their gravitational field (metric) power spectra which are in general functions of 
time. Here, we introduce the notations and equations that will be used in subsequent calculations. 

As before (see [3|, [5|, and references therein), we are working with perturbed Priedmann- 
Lemaitre- Robertson- Walker universes 

ds"^ = -c^dt^ + a^{t){5ij + hij)dx^dx^ = a^{'q)[~d-rf + {5ij + hij)dx'' dx^\, 

where the functions hij{ri,x.) are metric perturbation fields. Their spatial Fourier expansions are 
given by 

C /•+°° d^ 
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(1) 



The polarization tensors Pij (n) (s = 1,2) refer either to the two transverse-traceless components 
of gravitational waves (gw) or to the scalar and longitudinal-longitudinal components of density 
perturbations (dp). Density perturbations necessarily include perturbations of the accompanying 

matter fields (not shown here). 

s si" 

In the quantum version of the theory, the quantities c^ and Cn are the annihilation and creation 

operators, respectively, of the considered type of perturbations, and the |0) is the initial vacuum 

(ground) state of the corresponding time-dependent Hamiltonian. The metric power spectrum 



h'^{n,r]) is defined by the expectation value of the quadratic combination of the metric field: 

oo 
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The mode functions hn iv) ^r^ taken either from gw or dp equations, and C = \/16ttIp\ for 
gravitational waves and C = \/247r/pi for density perturbations. 

The simplest assumption about the initial stage of cosmological expansion (i.e. about the 
initial kick that presumably took place soon after the birth of our Universe [a, [3] ) is that it can be 



m 



described by a single power-law scale factor |7H9|] 

airj)=lo\v\'^^, (3) 

where lo and (3 are constants, /3 < — 1. Then the generated primordial power spectra (primordial 
means the interval of the spectrum pertaining to wavelengths longer than the Hubble radius at the 
considered moment of time) have the universal power-law dependence, both for gw and dp: 

/i2(n)cxn2(/5+2)_ 

It is common to write these power spectra separately for gw and dp: 

h^{n) (gw) = St^n"*, h^{n) (dp) = Sfn"'-^ (4) 

.^n accordance with the theory of quantum-mechanical generation of cosmological perturbations 
THSl, the spectral indices are approximately equal n^ — 1 = n^ = 2(/3 -|- 2) and the amplitudes B^ 
and Bg are of the order of magnitude of the ratio Hi/Hpi, where Hi ~ c/Iq is the characteristic 
value of the Hubble parameter during the kick. 

If the initial stage of expansion is not assumed to be a pure power-law evolution ([3]), the spectral 
indices rit and Ug — 1 are not constants, but their wavenumber dependence is calculable from the 
time dependence of the scale factor a{r]) and its time derivatives [lO|]. In fact, as we shall argue 
below, the CMB data suggest that even at a span of 2 orders of magnitude in terms of wavelengths 
the spectral index rig is not the same. We discuss this issue in detail in Sec. IIVI 

In what follows, we use the numerical code CAMB 

power spectra adopted there: 

^ \ "i / ^ \ "s-i 



111 ] and related notations for gw and dp 



m) = A,^-) , PM = A.^-J , (5) 

where /cq = 0.002Mpc~^. Technically, the power spectrum Ps{k) refers to the curvature perturba- 
tion called TZ or C, but the amplitudes Bg and {Ag)^''^ are equal to each other up to a numerical 



coefficient of order 1. The constant dimensionless wavenumber n is related to the dimensionful k by 
k = n/{2l}{), where Ih = c/Hq is the present-day Hubble radius. The wavenumber n = uh = 47r 
marks the wavelength equal to Ih today. The CMB temperature anisotropy at the multipole I 



is mostly generated by metric perturbations with wavenumbers n ^ i (see [13] for details). Set- 
ting h = 0.704 we obtain i w (0.85 x 10^Mpc)A;, which is consistent with the numerical result 



(1.0 X 10'^Mpc)fc derived in [l3(] 



The CMB temperature and polarization anisotropies are usually characterized by the four angu- 
lar power spectra: Cf , Cf , Cf , and Cj as functions of the multipole i. The contribution of 
gravitational waves to these power spectra has been studied, both analytically and numerically, in a 
number of papers 
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181 ]. The derivation of today's CMB power spectra brings us to approximate 
formulas of the following structure [iTI] : 

Cr = rithHn,Vrec)[Fr{n)]\ 

CF = I ^Hn,rjrec)h'{n,i^rec) [Fl{n)Ff{n)] , (6) 

Cr = I ^{h')\n,rjrec) [FY{n)]\ where Y = E,B. 
In the above expressions, h'^{n,r]rec) and {h')'^{n,rirec) are the power spectra of the gravitational 
wave field and its first time-derivative, respectively. The spectra are taken at the recombination 
(decoupling) time r]rec- The functions F^{n) {X = T,E,B) take care of the radiative transfer of 
CMB photons in the presence of metric perturbations. As already mentioned, the power residing in 
the metric fluctuations at wavenumber n mostly translates into the CMB TT power at the angular 
scales i ^ n. Similar results hold for the CMB power spectra induced by density perturbations. 
The actually performed numerical calculations use equations more accurate than Eq. ([6|). They 
also include the effects of the reionization era. 

The CMB power spectra needed for the analysis of WMAP data are calculated in the framework 
of background cosmological ACDM model characterized by the WMAP7 best-fit parameters [l|] 

fifo/i^ = 0.02260, r^c/i^ = 0.1123, Qa = 0.728, Trdon = 0.087, h = 0.704. (7) 

To quantify the contribution of relic gravitational waves to the CMB we use the quadrupole 
ratio R defined by 

CS(gw) 

Another measure is the so-called tensor-to-scalar ratio r. This quantity is constructed from the 
primordial power spectra (^ 

r = f. (9) 



Often this parameter is linked to incorrect (inflationary) statements, such as the inflationary 'con- 
sistency relation' (more details in Sec. lIVp . However, if one uses r without implying inflationary 
claims, one can find a useful relation between R and r. 

In general, the relation between R and r depends on background cosmological parameters and 
spectral indices n^ and n^. We found this relation numerically using the CAMB code [ll{]. (For a 
semianalytical approach see [19].) The results are plotted in Fig. [TJ For this calculation we used 
the background cosmological parameters ([7]) and the condition rit = Ug — 1 required by the theory 
of quantum-mechanical generation of cosmological perturbations. We verified that the relation 
r{R) only weakly depends on the background parameters and does not change significantly when 
the values d?]) are varied within the WMAP7 la error range [l(]. It is seen from the graph that 
r = 1.92R in the case of n^ = 1.0, and r ~ 2R for all considered n^, if R and r are sufficiently 
small. In other words, one can use r ~ 2R for a quite wide class of models. 

We do not know enough about the very early Universe to predict R with any certainty. However, 
since the theory of quantum-mechanical generation of cosmological perturbations requires that the 
amplitudes Bt and Bg (as well as At and As in Eq.([5])) should be of the same order of magnitude, our 
educated guess is that R should lie somewhere in the range R G [0.01, 1]. If i? were observationally 
found significantly outside this range, we would have to conclude that the underlying perturbations 
are unlikely to be of quantum-mechanical origin. On the other hand, the most advanced inflationary 
theories predict the ridiculously small amounts of gravitational waves, something at the level of 
r ~ 10~^^ or less, r G [0, 10^^^]. The rapidly improving CMB data will soon allow one to decisively 
discriminate between these theoretical frameworks (let alone the already performed discrimination 
on the grounds of purely theoretical consistency). We discuss this issue in Sec. IVIi 

III. EVALUATION OF RELIC GRAVITATIONAL WAVES FROM 7- YEAR WMAP DATA 

A. Likelihood function 

Relic gravitational waves compete with density perturbations in generating CMB temperature 
and polarization anisotropics at relatively low multipoles. For this reason we focus on the WMAP7 
data at 2 < ^ < 100. As before [sl, |]J, the quantities DJ^ , DJ^, Df^, and Df^ denote the 
estimators (and also the actual observed data points in the likelihood analysis) of the corresponding 
power spectra. Since the WMAP7 EE and BB observations are not particularly informative, we 
marginalize (integrate) the total probability density function (pdf) over the variables Df^ and 
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FIG. 1: The relation between R and r for different values of the spectral index Ug. From top to bottom the 
ns changes from Us = 0.8 to Us = 1.2. 
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y, bj. The resulting pdf [3| is a function of DJ'^ and DJ^: 

f{Dr,Dr)=n'x^-^{2'^-.T\'^){l-pl){ajr{aff] 



X exp 



P(.z 



2x{afY 2{ajy 



(10) 



This pdf contains the variables Df {XX' = TT, TE) through the quantities x = n{Dj + 
NJ ) and z = nDj , where NJ^ and Nf are total noise power spectra. Information about the 
power spectra Cf^ is contained in the quantities aj, af and pe (see [3| for details). The sought 
after parameters As, At, Ug, and nt enter the pdf through the Cf-^ . The quantity n = (2£ + l)/sky 
in (jlOp is the effective number of degrees of freedom at multipole i, where /sky is the sky-cut factor. 

In the WMAP7 data .lease .he sk,-o„t facto, . /.,, = 0.783 [2], wh.ch .s shghtly s.ahet than 
/sky = 0.85 used in WMAP5 data analysis '3"]. The smaller /sky increases the uncertainties, but 
this disadvantage is more than compensated by the reduction of overall noises. Therefore, the error 
bars surrounding the WMAP7 data points are somewhat smaller than those for the WMAP5 data 
release. This fact, together with slightly shifted data points themselves, allows us to strengthen 
our conclusions (see below) about the presence of gravitational wave signal in the WMAP data. 



We seek the perturbation parameters R, Ag and n^ (n^ = n^ — 1) along the hnes of our maximum 
Hkehhood analysis of WMAP5 data [3|]. The pdf (|10|) considered as a function of unknown R, As, 
and rig with known data points D^^ is a likelihood function subject to maximization. For a set 
of observed multipoles £ = 2, • • -jlmax, the likelihood function can be written as [3| 

where the constant C is chosen to make the maximum value of C equal to 1. 



B. Results of the analysis of the WMAP7 data 

The WMAP7 data points for DJ^ and DJ^ at multipoles 2 < i < l^nax = 100 were taken, 
with gratitude, from the Website [20[. The 3-dimensional likelihood function (jlip was probed by 



the Markov chain Monte Carlo method 



21 



22l | using 10,000 samples. The ML values of the three 



perturbation parameters R, Us and Ag were found to be 

R = 0.264, Us = 1.111, As = 1.832 x 10"^ (12) 

and nt = 0.111. In Fig. [2]we show the projection of the 10,000 sample points on the 2-dimensional 
planes R — Us and R — As. The color of an individual point signifies the value of the 3-dimensional 
likelihood of the corresponding sample. The projections of the maximum (I12p are shown by a black 
+. (The value R = 0.264 is equivalent to r = 0.550.) 

Before analyzing the 3-parameter results, it is instructive to consider the 2-parameter and 
1-parameter probability distributions of the sought after parameters. These marginalized distri- 
butions are obtained by integrating the likelihood function C pT]) (already represented by 10,000 
points) over one or two parameters. By integrating over As or n^, we arrive at 2-dimensional 
distributions for the pairs R — Ug ox R — As, respectively. The area around the maximum of the 
resulting distributions is shown in Fig. [3l In the R — Us space, the maximum is located at 

R = 0.228, Us = 1.108. (13) 

In the R — As space, the maximum is located at 

R = 0.253, As = 1.787 x 10~^. (14) 

In the left panel we also reproduce the 2-dimensional contours obtained by the WMAP team [l|, 2 1 
(their r is translated into our R). In contrast to the WMAP5 paper [23] . the WMAP7 paper 2| 







1 0.2 0.3 0.4 05 0.6 0.7 

Likelihood 



FIG. 2: The projection of 10,000 samples of the S-dimensional hkelihood function, based on the WMAP7 
data, onto the R — Us (left panel) and R — As (right panel) planes. The black + indicates the maximum 
likelihood parameters listed in (fT2|) . 
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FIG. 3: The ML point (red x) and the 68.3% and 95.4% confidence intervals (red solid lines) for the 
2-dimensional distributions R — Ug (left panel) and R — As (right panel). The left panel shows also the 
WMAP7 confidence contours (black dashed line) derived under the assumption that the spectral index Ug 
is one and the same constant throughout all measured multipoles [l|, l2|- The WMAP7 papers do not show 
the confidence contours in the R — As plane. 



shows only the uncertainty contours derived under the assumption of a strictly constant spectral 
index n^, but not for the case of its running. 

The 1-dimensionaI distributions for R, Ug or Ag are obtained by integrating the likelihood 
function £ (fTT|) over the sets of two parameters {As, n^), {As, R) or {R, Us), respectively. These 
distributions are presented in Figs. [Hand [S] (red solid lines). The ML values of the parameters 
and their 68.3% confidence intervals are found to be 

R = 0.273tm ns = l.ll2t',Z, ^ = (1.7651°;^^^) x lO^^. (15) 

For completeness and comparison (see Sec. lIIlCj) . we also show in Figs. [3] and [S]the 1-dimensional 
(1-D, for brevity) distributions derived by exactly the same procedure from WMAP5 and WMAP3 
observations. Obviously, the WMAP5 curves are copies of the previously reported distributions 
[31]; the WMAP3 curves are explained in Sec. IIIICl 

The constancy of spectral indices is a matter of special discussion in Sec. llVi In preparation 
for this discussion, we report the results of our likelihood analysis of data in other intervals of 
multipole moments i, in addition to the interval 2 < i < 100 that resulted in the ML values (|12p . 
First, we analyzed the WMAP7 data in the interval 101 < i < 220. The Ug coordinate of the 
maximum in 3-dimensional space R, Ag, and Ug was found to be ng = 0.951. The 1-D marginalized 
distribution for Ug gave the ML result 

ng = 0.969+|J:JJ| (68.3% C.L.). (16) 

Second, we have done the same in the combined interval of multipoles 2 < I < 220. The maximum 
of 3-D likelihood lies at n, = 1.003, whereas the 1-D ML result is 



Ug = 1.0211°:°^^ (gg 3% c L )^ (-17) 



The 1-D results allow one to make easier comparison of confidence intervals. As is seen from (J15p 
and (|16|) . the 1-D determinations of Ug in the adjacent ranges 2 < i < 100 and 101 < i < 220 
overlap, but only marginally, in their la intervals. As expected, if Ug is assumed constant in the 
entire range 2 < i < 220, its value (J17p is intermediate between (J15p and (J16p . The spreads of 
wavenumbers to which the 3-d spectral indices ng = 1.111, Ug = 0.951 and Ug = 1.003 refer are 
shown by marked red lines in Fig. [6l 

C. Comparison of results derived from WMAP7, WMAP5 and WMAP3 data releases 

One and the same 3-parameter analysis of WMAP5 and WMAP7 data has resulted in somewhat 
different ML parameters. From WMAP7 observations we extracted the ML parameters ()12p . 
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FIG. 4: The 1-dimensional likelihoods for R. The resuhs of the analysis of WMAP7, WMAP5 and WMAP3 
data are shown, respectively, by the red (solid), black (dashed), and blue (dotted) curves. The shaded 
regions indicate the 68.3% and 95.4% confidence intervals for WMAP7 likelihood. 

whereas the ML parameters extracted from WMAP5 observations [3] are 

R = 0.229, Us = 1.086, As = 1.920 x 10"^ 

Certainly, the results are consistent and close to each other. The same holds true for marginalized 
1-dimensional parameters and distributions shown in Figs. H] and [SI This similarity of results 
testifies to the stability of data and data analysis. There exist, however, important trends in the 
sequence of ML parameters extracted from the progressively improving WMAP3, WMAP5, and 
WMAP7 data. We want to discuss these trends. 

Specially for this discussion, we derived the parameters R,ns, and Ag from WMAP3 data in 
exactly the same manner as was done here and in [3] with WMAP7 and WMAP5 data releases. 



Previously 



12l |. we derived these parameters by a different method: we restricted the likelihood 



analysis to TE data and a single parameter R, while Ug and Ag were determined from phenomeno- 
logical relations designed to fit the TT data. That analysis has led us to i? = 0.149]'^Q']^4g and 
ng = 1.002. Our new derivation, based on 3-dimensional likelihood, gives the following WMAP3 
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FIG. 5: The 1-dimensional likelihoods for Ug (left) and Ag (right). In both panels the red (solid), black 
(dashed) and blue (dotted) curves show the results for WMAP7, WMAP5 and WMAP3 data, respectively. 

maximum likelihood parameters: 

i? = 0.181, n, = 1.045, ^, = 2.021 X 10"^ 



The corresponding 1-dimensional distributions give 



R = 0.205 



+0.181 
-0.157' 



n, = 1.059 



h0.097 
-0.066' 



As = {l.SMtlil'l) X 10"^ 



(18) 



The 1-dimensional WMAP3 distributions are plotted in Figs. [Hand [5] by dotted curves alongside 
with WMAP5 (dashed) and WMAP7 (solid) curves. 

Looking at all derivations and graphs collectively, we can draw the following conclusions. First, 
about R. The maximum likelihood value of R increases as one goes over from 3-year to better 
quality 7-year data. The ML value R = 0.264 of WMAP7 data release is 15% larger than the 
analogous result R = 0.229 obtained from WMAP5 data. The 3-, 2-, and 1-parameter determina- 
tions of R persistently concentrate somewhere around the mark R = 0.25. The uncertainties AR, 
although still considerable, get smaller as one progresses from 3-year to better quality 7-year data. 
The R = (no gravitational waves) hypothesis is under increasing pressure. For example, the 
WMAP7 1-parameter result, Eq. (fT5]) and Fig. [H excludes the R = hypothesis at almost 2a level. 
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To be more precise, the R = point is right on the boundary of the 94% confidence area (94% of 
surface area under the WMAP7 curve in Fig. HD surrounding the 1-D ML value R = 0.273. This 
is an improvement in comparison with a shghtly larger than Icr interval separating R = from 
the WMAP5 1-d ML value R = 0.266 [3]. Admittedly, the gradual changes in ML values of R are 
small, while uncertainties are still large. It would not be very surprising if random realizations of 
noise moved the consecutive 3-year, 5-year, 7-year ML values of R in arbitrary order. Nevertheless, 
we observe the tendency for systematic increase and saturation of R alongside with decrease of 
AR. 

Second, our conclusions about Ug. Together with the increase of ML R, we observe the tendency 
for the increase of n^ accompanied by the theoretically expected and understood decrease of Ag 
(the contribution of relic gravitational waves becomes larger, while the contribution of density 
perturbations becomes smaller). The ML value Ug = 1.111 derived from the WMAP7 data, Eq. p2p . 
is larger than Ug = 1.086 derived from WMAP5 data, while the WMAP7 value of Ag is somewhat 



Is 



smaller than that of WMAP5 [31]. The tendency for increase of Ug and decrease of Ag is also 
illustrated by the 1-parameter distributions in Fig. [5l The spectral indices ng,nt persistently point 
out to blue primordial spectra, i.e. Ug > 1 for density perturbations and nj > for gravitational 
waves, in the interval of wavelengths responsible for CMB anisotropics at 2 < ^ < 100. The 
larger values of Ug, nt derived from WMAP7 data release enhance the doubt on whether the 
conventional scalar fields could be the driver of the initial kick, since these cannot support /3 > —2 
in Eq.([3|) and, consequently, n^ > 1, nj > in Eq.([l]). (Since the inflationary theory is capable of 
predicting virtually anything that one can possibly imagine, there exist of course literature claiming 
that inflation can predict blue spectra. But these claims are based on the incorrect (inflationary) 
formula for density perturbations, see [9|.) The questions pertaining to the spectral indices are 
analyzed in some detail in the next section. 

IV. COMPARISON OF OUR RESULTS WITH CONCLUSIONS OF WMAP COLLABO- 
RATION 

Having analyzed the 7-year data release, the WMAP collaboration concludes that a minimal 
cosmological model without gravitational waves and with a constant spectral index Ug across the 
entire interval of considered wavelengths remains a "remarkably good fit" to ever improving CMB 
data and other datasets. The WMAP team emphasizes: "We do not detect gravitational waves 
from inflation with 7-year WMAP data, however the upper limits are 16% lower..." [2], P-H; "The 
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7-year WMAP data combined with BAO and Hq excludes the scale-invariant spectrum by more 

than 3o", if we ignore tensor modes (gravitational waves)" [l], p. 15; "We find no evidence for tensor 

i i fl 

modes..." [1!], p. 30; "We find no convincing deviations from the minimal model" [ll], p.l, etc. 

In contrast, our analysis of WMAP3, WMAP5 and WMAP7 data leads us in the opposite 
direction: the improving data make the gw indications stronger. The major points of tension 
between the two approaches seem to be the constancy of spectral indices and the continuing use 
by the WMAP team of the inflationary theory in data analysis and interpretation. We shall start 
from the discussion of spectral indices. 

The constancy of spectral indices is a reasonable assumption, but not a rule. If the power-law 
dependence ([3]) is not a good approximation to the gravitational pump field during some interval 
of time, the constancy of rit, Us is not a good approximation to the generated primordial spectra 
^ in the corresponding interval of wavelengths. In fact, the future measurements of frequency 
dependence of the spectrum of relic gravitational waves will probably be the best way to infer the 
"early history of the Hubble parameter" [lOJ] . 

The frequency-dependence of a gw spectrum is fully determined by the time-dependence of 
the function ^{t) = —H/H^. (In more recent papers of other authors this function is often 
denoted e(t).) The function j{t) describes the rate of change of the time-dependent Hubble radius 
Init) = c/H{t): 

(f) = ± (J—\ = \dhi(^ 
^^ ' dt \H{t)) c dt • 

The function 7(t) is a constant for power-law scale factors ^: 7 = (2 + /3)/(l + /3), and 7 = for 
a period of de Sitter expansion. The interval of time dt during the early era when gravitational 
waves were entering the amplifying regime and their today's frequency spread dv are related by 
(see Eq.(21) in [10(): 



dt d\r\.i/ 

The today's spectral energy density of gravitational waves e{v) is related to the early universe 

parameter 7(t) by (see Eq.(22) in [lOl]): 

[d \iie{v)/d \nv] 
^^"^ = -2-[dlne{.)/dln.y ^^^^ 

The spectral index Ug of a pure power-law energy density e(i^) oc u"^ is defined as Ug = 
[d In e{iy)/d Inz^]. It is reasonable to retain this definition for more complicated spectra. Then, 
Eq. (119p can be rewritten as 

.»(.) ^ -^. (20, 
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Obviously, in the case of pure power-laws ([3]) we return to the constant spectral index Ug = 
-27/(1 -7) = 2(2 + /3). 

Equation ()20p was derived for the energy density of relatively high-frequency gravitational 
waves, i^ > 10^^^ Hz, which started the adiabatic regime during the radiation-dominated era. In 
our CMB study we deal with significantly lower frequencies. It is more appropriate to speak 
about wavenumbers n rather than frequencies i^, and about power spectra of metric perturbations 
h'^{n) rather than energy density e(z^). The ^-dependent spectral index nt{k) entering primordial 
spectrum ([5]) is defined as nt{k) = [d In Pt{k)/d Ink]. Then formula for nt{k) retains exactly the 
same appearance as Eq. (pO|) : 

The spectral index nt{k) reduces to a constant rit = 2(2 -|- /3) in the case of power-law functions 

The spectral index ns{k) — 1 for density perturbations is defined by ns{k) — 1 = 
[d lnPs{k)/d Ink]. The formula for ns{k) — 1 is more complicated than Eq. (j2ip as it contains 
also dj{t)/dt as a function of k. However, it is important to stress that adjacent intervals of 
power-law evolution ([3]) with slightly different constants /3, will result in slightly different pairs of 
constant indices rit, n^ — 1 in the corresponding adjacent intervals of wavelengths. Of course, the 
spectrum itself is continues at the wavelength marking the transition between the two regions. 

Extending the minimal model, the WMAP Collaboration works with the power spectrum [l| 

, \ n3(fco)-l+5«3ln(fc/fco) 



Ps{k) = Psiko) ( - 1 , (22) 

which means that the /c-dependent (running) spectral index ns{k) is assumed to be a constant plus 
a logarithmic correction: 

ns{k) = ns{ko) + as ln{k/ko). (23) 

The aim of WMAP data analysis is to find a^, unless it is postulated from the very beginning, 
as is done in the central (minimal) model, that as = 0. We note that although logarithmic 
corrections do arise in simple situations and can even be termed "natural" [lO(], they are not 
unique or compulsory, as we illustrated by exact formula (j2ip . Nevertheless, we do not debate this 
point. We accept WMAP's definitions, and we want to illustrate their results graphically, together 
with our evaluations of n^ in this paper. 

The main result of WMAP7 determination is Ug = 0.963 it 0.012 (68% C.L.) derived under 
the assumption of no gravitational waves and constant n^ throughout all wavelengths included 
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in the considered datasets [1|, y]. When the presence of gravitational waves is allowed, but n^ is 
still assumed constant, the n^ rises to n^ = 0.982^Q'Q-^g from WMAP7 data alone. Finally, from 
WMAP7 data alone, the WMAP team finds ns{ko) = 1.027, a^ = -0.034 in the case of no gw 
but with running of n^, and ns(A;o) = 1.076, a^ = —0.048 in the case of running and allowed 
gravitational waves (we quote only central values without error bars, see Table 7 in [1]). All the 
resulting values of ns{k) derived by WMAP team are shown in Fig. [6j For comparison, we also 
plot by red lines our evaluations of n^, see Sec. IIIIBI 

The lines in Fig. [6] show clearly that our finding of a blue shape of the spectrum, i.e. Ug = 1.111, 
at longest accessible wavelengths is pretty much in the territory of WMAP findings, if one allows 
running, even as simple as Eq. (|23p . and especially when running is combined with gravitational 
waves. On the other hand, as was already explained in [^, the attempt of constraining relic 
gravitational waves by using the data from a huge interval of wavelengths and assuming a constant 
Ug (or its simple running) across all wavelengths is unwarranted. The high-£ CMB data, as well as 
other datasets at relatively short wavelengths, have nothing to do with relic gravitational waves, and 
their use is dangerous. As we argued in Sec. IIIIBt the spectral index Ug appears to be sufficiently 
different even at the span of two adjacent intervals of wavenumbers. The restriction to a relatively 
small number of multipoles 2 < £ < 100 is accompanied by relatively large uncertainties in R, but 
there is nothing we can do about it to improve the situation, this is in the nature of efforts aimed 
at measuring R. The difference in the treatment of Ug is probably the main reason why we do see 
indications of relic gravitational waves in the data, whereas the WMAP team does not. 

Another contributing factor to the difference of conclusions is the continuing use by WMAP 
Collaboration of the inflationary theory and its (incorrect) relation nt = —r/8, which automatically 
sends r to zero when rit approaches zero. This formula is a part of the inflationary 'consistency 
relations' 

r = 16e = —Snt- 

Only one equality in this formula, 16e = —8nt, is correct being an approximate version (for small 
7, e = 7) of our exact formula, Eq. (|2ip . The 'consistency relation' r = 16e is incorrect. It is 
an immediate consequence of the "classic result" of inflationary theory, namely, the prediction of 
arbitrarily large amplitudes of density perturbations generated in the limit of de Sitter inflation 
(e = 0, Ug = 1), regardless of the strength of the generating gravitational field (curvature of 



space-time) regulated by the Hubble parameter H 



26(1 . Certainly, it would be inconsistent, even 



by the standards of inflationary theory, not to use the relation r = —8nt in data analysis, if the 
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FIG. 6: The spectral index Us as a function of the wavenumber k. The ML results of this work are shown 
by red lines. Other lines are our plots of WMAP7 findings [l|, Table 7. 

inflationary "classics" is used in derivation of power spectra and interpretation of results (see, for 
example, Fig. 19 in [l]). Obviously, in our analysis we do not use the inflationary theory and its 
relations. 



V. FORECASTS FOR THE PLANCK MISSION BASED ON THE RESULTS OF ANAL- 
YSIS OF WMAP7 DATA 

The Planck satellite [4| is currently making CMB measurements and is expected to provide data 
of better quality than WMAP. We hope that the indications of relic gravitational waves that we 
found in WMAP3, WMAP5 and WMAP7 data will become a certainty after Planck observations. 
We quantify the detection ability of the Planck experiment by exploring the vicinity of the WMAP7 
maximum likelihood parameters (I12p . 

It is seen from Fig. [2] that the samples with relatively large values of the likelihood (red, yellow, 
and green) are concentrated along the curve which projects into relatively straight lines (at least, 
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up to i? ~ 0.5) in the planes R — Ug and R — Ag'. 

Hs = 0.98 + 0A9R, As = (2.30 - 1.77i?) x 10"^ (24) 



The ML model ()12p is a specific point on these lines, corresponding to i? = 0.264. The pa- 
rameterization (j24p is close to the result derived from WMAP5 data: Ug = 0.98 + 0.46i?, Ag = 
(2.27 - 1.53R) X 10^9, and ML R = 0.229 (see Eq. (15) in [3]). We use Eq.dMD in formulating our 
forecast, thus reducing the task of forecasting to a 1-parameter problem in terms of R. 
Following [31] , we define the signal-to-noise ratio as 

S/N = R/AR, (25) 

where the numerator is the true value of the parameter R (or its ML value, or the input value in a 
numerical simulation) while AR in the denominator is the uncertainty in determination of R from 
the data. 

We estimate the uncertainty AR using the Fisher matrix formalism. We take into account all 
available information channels, i.e. TT, TE, EE and BB correlation functions, and their various 
combinations. The uncertainty AR depends on instrumental and environmental noises, on the 
statistical uncertainty of the CMB signal itself, and on whether other parameters, in addition to R, 
are derived from the same dataset. All our input assumptions about Planck's instrumental noises, 
number and specification of frequency channels, foreground models and residual contamination, 
sky coverage and lifetime of the mission, etc. are exactly the same as in our previous paper [3]. We 



do not repeat the details here and refer the reader to the text and appendices in 



3] which contain 



necessary references. Technically, our present forecast is somewhat different (better) than that in 
the previous analysis [31] because of slightly different family of preferred perturbation parameters 



(j24p . with slightly higher than before the ML value of i?, R = 0.264. We have added only one 
new calculation, reported in Fig. [11] and Fig. W2\ which is the S/N for the survey of 28 months 
duration, instead of the nominal assumption of 14 months. 

The results of our forecast for the Planck mission are presented in figures. In Fig. [7] we show 
the total S/N for the case TT + TE + EE + BB, i.e. when all correlation functions are taken into 
account, and at all relevant multipoles 2 < i < 100. The possible levels of foreground cleaning 
are marked by a^^. The pessimistic case is the case of no foreground removal, a^^ = 1, and the 
nominal instrumental noise of the BB information channel at each frequency increased by a factor 
of 4. Three frequency channels at lOOGHz, 143GHz and 217GHz are considered as providing data 
on perturbation parameters -R, n^, Ag. The more severe. Dust A, model is adopted for evaluation 
of residual foreground contamination. 
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FIG. 7: The total signal-to-noise ratio SjN as a function of i?. The left panel shows SjN with Ai? 



^/VPrb.: while the right panel shows 5/A^ with Ai? = \/iF ^)b.r- 



In the left panel of Fig. [7] we consider the idealized situation, where only one parameter R 
is unknown and therefore the uncertainty Ai? is calculated from the Ffm element of the Fisher 



matrix, Ai? = 1/^/Frr. In the right panel of Fig. [7] we consider a more realistic situation, where 
all perturbation parameters R, Us, Ag are unknown and therefore the uncertainty Ai? increases 
and is calculated from the element of the inverse matrix, Ai? = y^(i^~^)/j/j. 

From the right panel of Fig. [7] follows our main conclusion: the relic gravitational waves of the 
maximum likelihood model (J12p will be detected by Planck at the healthy level S/N = 4.04, even 
in the pessimistic case. This is an anticipated improvement in comparison with our evaluation 
S/N = 3.65 based on WMAP5 data analysis pj. The detection will be more confident, at the 
level S/N = 7.62, 6.91, if a^^ = 0.01, 0.1 can be achieved. Even in the pessimistic case, the 
signal-to- noise ratio remains at the level S/N > 2 for R > 0.11. 

Further insight in the detection ability of Planck and interpretation of future results is gained 
by breaking up the total S/N into contributions from different information channels and individual 
multipoles. It is easier to do this for the idealized situation, Ai? = 1/\/Fr^, exhibited in the left 
panel of Fig. [71 In Fig. [5] we show how the TT + TE + EE and BB contribute to the total S/N 
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FIG. 8: The decomposition of the total SjN into contributions from different information channels. Four 
panels describe different assumptions about foreground cleaning and noises. The SjN line from the upper 
left panel for the optimistic case a"^^ = 0.01 is copied as a broken line in other panels. 

based on all correlation functions TT + TE + EE + BB. (The (S/N)'^ for the full combination 
TT + TE + EE + BB is the sum of {S/Nf for TT + TE + EE and BB.) It is seen from the upper 
left panel of Fig. [8] that the i?-mode of polarization is a dominant contributor to the total S/N only 
in the conditions of very deep cleaning, o"^ = 0.01, and relatively small values of the parameter 
R. On the other hand, in the pessimistic case, the BB channel provides only S/N = 2.02 for the 
benchmark case R = 0.264. Most of the total S/N = 7.32 in this case comes from the TT+TE+EE 
combination. 

In Fig. [9] we illustrate the decomposition of the total S/N into contributions from individual 
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FIG. 9: The individual terms {S/N)'j as functions of ^ for various combinations of information channels and 
two opposite assumptions about residual foreground contamination and noises. The calculations are done 
for the ML model ([121) with R = 0.264. 



multipoles t. 



{S/Nf = Y,{S/N)l 



(26) 



We show the contributing terms {S/N)'j for three combinations of information channels, TT + 
TE + EE + BB, TT + TE + EE and BB alone, and for two opposite extreme assumptions about 
foreground cleaning and noises, namely, a^ = 0.01 and the pessimistic case. The calculations are 
done for the benchmark model (I12p with R = 0.264. Surely, the total S/N exhibited in the upper 
left and lower right panels of Fig. [8] is recovered with the help of Eq. (l26p from the sum of the terms 
{S/N)'j shown in the left and right panels of Fig. O respectively. 

It is seen from the left panel of Fig. [U] that in the case of deep cleaning the BB channel is 
particularly sensitive to the very low multipoles £ ~ 10 associated with the reionization era. At the 
same time, the right panel of Fig. [9] demonstrates that in the pessimistic case most of S/N comes 
from TT + TE + EE combination and, specifically, from the interval of mutlipoles i ~ (20 — 60) 
associated with the recombination era. 

Finally, we want to discuss possible improvements in our forecasts. They will be achievable, if 7 
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FIG. 10: The improved signal-to-noise ratio SjN for the scenario where the Dust B model is correct and 
the instrumental noises are smaller because of the 7 frequency channels used, instead of 3, in data analysis. 



The left panel shows SjN with Ai? == 1/\/Frr, while the right panel shows S/N with Ai?^ = ^/{F^^rb.. 



frequency channels, instead of 3, can be used for data analysis, and/or if the less restrictive Dust B 
model, instead of the Dust A model, turns out to be correct, and/or if 28 months of observations, 
instead of 14 months, can be reached. In Fig. [10] we show the results for the scenario where 7 
frequency channels are used and the Dust B model is correct. Similarly to Fig. [TJ the left panel 
shows S/N calculated assuming that only R is being determined from the data, whereas the right 
panel shows a more realistic case in which all perturbation parameters are unknown. In the later 
case, the S/N for the ML model ()12p with R = 0.264 increases up to the level S/N = 5.56 as 
compared with S/N = 4.04 that we found in the right panel of Fig. [7| for 3 frequency channels and 
the Dust A model. 

The two panels in Fig. [11] illustrate the improvements, as compared with the two panels in 
Fig. [T] arising from the longer, 28-month, survey instead of the nominal 14 months survey. The 
S/N for the benchmark model R = 0.264 increases up to S/N = 5.39 even in the pessimistic case, 
right panel. The decomposition of S/N appearing in the left panel of Fig. [11] into contributions 
from different information channels is presented in Fig. [12] One can see that the relative role 
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FIG. 11: The improved SjN for 28 months of observations. Other assmiiptions are the same as in Fig. [71 



of the BB channel has increased. (The prospects of i?-mode detection in the conditions of low 



instrumental and foreground noises have been analyzed in 



24|. 



VI. THEORETICAL FRAMEWORKS THAT WILL BE ACCEPTED OR REJECTED BY 
PLANCK OBSERVATIONS 

Our forecast for the Planck mission, as any forecast in nature, can prove its value only after 
actual observation. We predict sunny days of confident detection of relic gravitational waves, but 
the reality can turn out to be gloomy days of continuing uncertainty. One can illustrate this 
point with the help of probability distributions in Fig. HI We expect that Planck observations will 
continue the trend of exhibiting narrower likelihoods with maximum in the area near the WMAP7 
ML value R = 0.264. Then the detailed analysis in Sec. IVl explains the reasons for our optimism. 
But, in principle, the reality can happen to be totally different. From the position of pure logic it 
is still possible that the Planck data, although making the new likelihood curve much narrower, 
will also shift the maximum of this curve to the point i? = 0. In this case, instead of confident 
detection we will have to speak about "tight upper limits" . We do not think, though, that this is 
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FIG. 12: Decomposition of the total SjN figuring in the left panel of Fig. [TT] into contributions from 
different information channels. In the assumed conditions of deep cleaning a^^ = 0.01, the BB alone is more 
informative than the combination TT + TE + EE. 

going to happen. What kind of conclusions about theoretical models can we make if the likelihood 
curve conies out as we anticipated? 

As was already mentioned in Sec. m the theory of quantum-mechanical generation of cosmo- 
logical perturbations implies a reasonable guess for the true value of i?: R £ [0.01, 1]. We shall 
call it a model Mi. At the same time, the most advanced string-inspired inflationary theories 
predict R somewhere in the interval R G [0, 10^^^]. We shall call it a model M2. The inflationary 
calculations can be perfectly alright in their stringy part, but the observational predictions are 
entirely hanging on the inflationary "classic result", and therefore should fall together with it, on 
purely theoretical grounds. But this is not the point of our present discussion. We wish to conduct 
a Bayesian comparison of models Mi and M2, regardless of the motivations that stayed behind 
these models. 

The model Mi suggests that the quadrupole ratio R should lie in the range [0.01, 1] with a 
uniform prior in this range: Pprior(-R|Mi) = ci for R G [0.01,1]. The model M2 suggests that R 
should lie in the range [0, 10~^^] with a uniform prior in this range: Pprior(-R|M2) = C2 for R G 
[0, 10~^^]. The constants ci and C2 are determined from normalization of the prior distributions, 
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ci = (0.99)~^ and C2 = 10^^. The predicted interval of possible values of R is much wider for Mi 
than for M2. Therefore, Mi is penalized by a much smaller normalization constant ci than 02- 
The observed data allow one to compare the two models quantitatively with the help of the Bayes 
factor Ku uA 

^ f,^ PpnoriR\M,)C{R)dR 
'' J^ P^nor{R\M2)C{R)dR' ^ > 

where C is the likelihood function of R derived from the observation. 

We shall start from the existing WMAP7 data. The likelihood function as a function of R is 
shown by a red solid line in Fig. HI Calculating the Bayes factor according to Eq. (j27p we find 
K12 = 1.61, i.e. Ini^i2 = 0.48. As expected, this value of the Bayes factor, although indicative, 
does not provide sufficient reason for the rejection of model M2 in favor of Mi. 

With more accurate Planck observations the situation will change dramatically, if the WMAP7 
maximum likelihood set of parameters (J12p is correct. To make calculation in Eq. (j27p . we adopt a 
Gaussian shape of the likelihood function C with maximum at R = 0.264 and i? > 0. The standard 
deviation is taken as Ai? = 0.065. This value follows from the analysis of S/N for Planck mission 
and corresponds to the derived S/N = 4.04 in the pessimistic case. Then the calculation of K\2 



m, 



gives the value i^i2 = 579.23, i.e. Ini^i2 = 6.36. In accordance with Jeffrey's interpretation 
this result shows that the Planck observations will decisively reject the model M2 in favor of Mi. 
If Planck observations are as accurate as expected, and if our assumptions about Planck's 
likelihood function are correct, the models much less extreme than M2 will also be decisively 
rejected. We introduce the model M3 with a flat prior in the range R G [0, x] and ask the question 
for which x the Bayes factor iTis exceeds the critical value iiTia = 100 (Ini^Tia = 4.61) [25] which is 
required for decisive exclusion of the model. The calculation gives x = 0.05. This means that under 
the conditions listed above the Planck experiment will be able to decisively reject any theoretical 
framework that predicts i? < 0.05. In terms of the parameter r this means the exclusion of all 
models with r < 0.095. 

VII. CONCLUSIONS 

The analysis of the WMAP7 data release amplifies observational indications in favor of relic 
gravitational waves in the Universe. The WMAP3, WMAP5, and WMAP7 temperature and 
polarization data in the interval of multipoles 2 < i < 100 persistently point out to one and 
the same area in the space of perturbation parameters. It includes a considerable amount of 
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gravitational waves expressed in terms of the parameter R = 0.264, and somewhat blue primordial 
spectra with indices n^ = 1.111 and n^ = 0.111. If the maximum likelihood set of parameters that 
we derived from this analysis is a fair representation of the reality, the relic gravitational waves will 
be detected more confidently by Planck observations. Even under pessimistic assumption about 
hindering factors, the expected signal-to-noise ratio should be at the level S/N = 4.04. 
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